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Abstract 

"Mayer waves" are long-period (6 to 12 seconds) oscillations in ar- 
terial blood pressure, which have been observed and studied for more 
than 100 years in the cardiovascular system of humans and other mam- 
mals. A mathematical model of the human cardiovascular system is 
presented, incorporating parameters relevant to the onset of Mayer 
waves. The model is analyzed using methods of Lyapunov stabil- 
ity and Hopf bifurcation theory. The analysis shows that increase in 
the gain of the baroreflex feedback loop controlling venous volume 
may lead to the onset of oscillations, while changes in the other pa- 
rameters considered do not affect stability of the equilibrium state. 
The results agree with clinical observations of Mayer waves in human 
subjects, both in the period of the oscillations and in the observed 
age-dependence of Mayer waves. This leads to a proposed explana- 
tion of their occurrence, namely that Mayer waves are a "gain-induced 
instability" . 

1 INTRODUCTION 

The existence of fluctuations in blood pressure has been know since the intro- 
duction of the recording manometer by C. Ludwig, see [18]. These fluctua- 
tions, usually referred to as waves, are classified by various methods including 
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the name of the discoverer, the origin, physiological cause, or the frequency 
or period. The term Mayer waves refers to periodic fluctuations in blood 
pressure which are slower than respiration in animals with normal respira- 
tory movements. They were announced by S. Mayer in 1876 [15], and hence 
the name. They are also known as third order waves. 

The frequencies reported by various authors for Mayer waves differ con- 
siderably [18]. Those described by Mayer in rabbits had a frequency of 6- 
9 waves/min., while other researchers have found waves with frequencies 
ranging from 7-12 waves/min. in humans [18]. Some researchers have pro- 
posed to designate these waves as the "10-second-rhythm" [18]. The onset of 
Mayer waves may result in serious physiological implications, such as faint- 
ing. Mayer waves are of interest to researchers seeking to fully understand 
the functioning of the cardiovascular system. 

It is generally conceded that Mayer waves appear most often when the 
subjects are exposed to abnormal conditions [2]. Lack of oxygen, the effects 
of severe haemorrhage, and other extreme or sudden changes in blood supply 
to parts of the body favour the appearance of these slow periodic fluctuations 
[2]. Experiments have shown that when the blood pressure is measured for 
subjects lying in a supine position and then in a tilted position, there may 
exist Mayer-like oscillations for the tilted position. A remarkable feature 
observed in these experiments is that the Mayer waves occur more frequently 
in younger subjects, and disappear with age [6, 10, 12, 14, 16, 17]. The origin 
of Mayer waves, however, remains an unsolved problem [18]. At a conference 
on Mayer waves held in Prague in 1977, various theories were proposed to 
explain the origin of these waves. Four main theories, the myogenic theory, 
the central theory, the feedback theory, and the resonance theory were given. 

The myogenic theory states that third order waves are due to the inherent 
property of peripheral systems, (vascular smooth muscle), to exhibit sponta- 
neous rhythmic activity. The central theory postulates that cardiovascular 
centres in the brain stem generates the slow rhythms in a way similar to the 
respiratory centre. The feedback theory, as its name suggests, attributes the 
origin of Mayer waves to delays and nonlinearities in the body's feedback 
control mechanisms. The resonance theory on the other hand, postulates 
that one or more of the above factors enable the overall circulatory system 
to resonate at certain frequencies. Among the four, the feedback theory 
and the central theory are considered most probable [11, 21, 22]. However, 
the stability of the central oscillator is questionable [21]. Thus among the 
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two most popular theories, the feedback theory remains the most probable 
explanation for the existence of Mayer waves. 

It is essential, for survival, that blood pressure be controlled to stay within 
a narrow, safe range. This function is performed by the body's control mecha- 
nisms, the fastest being the baroreflex. The baroreceptors are stretch sensors 
located in the systemic arteries which detect changes in blood pressure. The 
baroreflex feedback loops respond to baroreceptor impulses to control blood 
pressure via three mechanisms: heart rate, systemic capillary resistance and 
venous volume. All three mechanisms are explored in this work. 

DcBoer et. al.[5] proposed that blood pressure variability is caused by a 
time delay in the baroreceptor loop. If this were the case, one would expect 
the delay to increase with age due to a slowing of the body's responses. 
This would then cause an increase in the existence of Mayer waves in older 
people. While we do not discount the contribution of delay, it appears not 
to be the main cause of Mayer waves. We hypothesize that blood pressure 
variability may be attributed to a change in feedback gain, extending the 
work of Wesseling et. al. [21] [22] . Previous studies of feedback-control systems 
in physiology (Glass and Mackey model) [7] [13], and in engineering (Watt s 
regulator model) [8, 19], have shown that an increase in feedback gain can 
cause a system to change behaviour from a steady state to an oscillating 
state. This may be called a "gain-induced instability" and has been studied 
by use of the Hopf Bifurcation Theorem. Since young adult humans tend to 
have quicker reflexes and better muscle tone than the elderly, they can be 
expected to have higher gain in the baroreceptor loop. Thus our hypothesis 
that Mayer waves may be a gain-induced instability is consistent with the 
observed age-related data. 

A mathematical model can be used to give greater insights into the roles 
of the various mechanisms affecting Mayer waves. Thus the primary objective 
of this study is to develop a dynamical model for the mammalian circula- 
tory system and use it to analyze blood pressure variability dynamics, as a 
function of the physiological parameters in the model. The dynamical model 
is a generalization of the steady-state model of Hoppensteadt and Peskin 
[9]. In addition to incorporating temporal dynamics, this model will allow 
investigation of the effects of each of the three baroreflex feedback loops, in- 
dependently of the other two. This type of experiment has been performed 
on animals, eg. dogs [20], but is difficult to carry out on human subjects. 
Parameter values in the mathematical model are chosen to correspond to a 
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typical adult human being. 

2 MODELLING 

While it is desirable to include the behaviour of each cardiovascular com- 
ponent in a model of the circulation, certain components can be lumped 
together without sacrificing the qualitative behaviour of the system [3, 14]. 
This section presents the modelling assumptions and the development of the 
model, first for the basic fluid flow of blood in the cardiovascular system, 
then with the nonlinear baroreflex control. 

2.1 MODELLING ASSUMPTIONS 

The assumptions and simplifications underlying the mathematical model are 
stated in the following. 

1. The cardiovascular system is a closed-looped hydrodynamic system 
comprising two heart pumps (the left and right sides of the heart), 
two large arteries, two veins and the two capillary networks, corre- 
sponding to the systemic and pulmonary circulations respectively. The 
total blood volume is constant in time. 

2. The large arteries and veins and the heart are comphance vessels [9], 
that is, volume is proportional to pressure in these vessels. On the other 
hand, the smaller arteries and veins in the capillary networks arc resis- 
tance vessels, that is, flow is proportional to pressure. The unstressed 
volume of blood vessels is negligible at all parts of the circulation except 
in the systemic veins. 

3. Flow from the heart is continuous, that is the pulsatile nature of blood 
pressure is neglected and only average pressures and volumes, over the 
period of the pulse, are dealt with. 

4. The pressure in the heart when relaxed, is equal to that of the veins 
supplying blood to it. That is, the pressure in the right and left hearts 
are those of the systemic and pulmonary veins respectively. On its 
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expansion stroke (diastole) the heart receives a volume of blood pro- 
portional to this venous pressure. The heart pumps out at each con- 
traction stroke (systole) the amount of blood received from the veins 
on the previous diastole. This is the Prank-Starhng model of the heart 



5. Cardiovascular blood pressure is controlled by a baroreceptor feedback 
mechanism, which acts on systemic venous unstressed volume, sys- 
temic venous compliance, systemic capillary resistance and heart rate, 
to counteract changes in systemic arterial pressure. 

6. The baroreceptor feedback gain and its dependence on systemic arterial 
pressure is modelled as a Hill function (described below). 

7. Changes in venous volume, systemic resistance and heart rate act in- 
dependently in parallel on blood pressure. 

8. Compliance is constant in all parts of the cardiovascular system except 
the systemic veins, where it may be varied by the baroreflex. 

9. Resistance is constant for the capillary networks of the pulmonary cir- 
culation, but may be varied by the baroreflex in the systemic circula- 



A notational convention adopted throughout this model is that dynamic 
variables are represented by lower case letters, while parameters and labels 
are written in upper case. We first construct a simple linear model of the 
cardiovascular system, then add the baroreflex control system. 

2.2.1 Linear Cardiovascular Model 

The following linear relationship between volume, t), pressure, p and compli- 
ance, C, in the large vessels (arteries and veins) of the circulation [9] , is the 
mathematical form of Assumption 2. (i.e. that these are compliance vessels): 



[22]. 



tion. 



2.2 



MODEL DEVELOPMENT 



V — C * p 



(1) 
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However, Equation (1) suggests that if p = 0, then v = 0, which is not the 
case. A more reahstic relation is: 



v^Vd + {C*p) (2) 

where Vd is the unstressed volume, that is the volume of the vessels ai p — 
0. The unstressed volume in the systemic venous circulation is very impor- 
tant, as over one-half of the venous volume is unstressed volume [4]. (see 
Assumption 2.). 

The flow, q in the vessels of the capillary networks, modelled as in [9] as 
resistance vessels, is (by Assumption 2): 

Pa - Pv .ON 
^^^^ 

Here pa, Pv and R represent the pressure in the arteries and veins and 

resistance respectively, in cither the systemic or pulmonary circulation. 

From the Frank-Starling Assumption 4, the following relations for the left 
and right cardiac outputs, and q'r, respectively are given: 

= F * Cl * Ppv ^ Kl* Ppv (4) 

qji^ F*Cr* psv ^ Kr* psv (5) 

Cl and Cr are the compliances in the left and right hearts respectively and 
F is the heart beat frequency. The subscripts A and V represent arteries 
and veins respectively, while the subscripts S and P stand for the systemic 
and pulmonary circulations, respectively. 

The rate of change of volume of an incompressible fluid in a vessel is the 
difference between the rates of flow of the fluid, into and out of the vessel. 
Hence, the following differential equations are obtained, for the change of 
volume of blood in the systemic arteries, systemic veins, pulmonary arteries 
and pulmonary veins respectively: 

dvsA /„v 
-^ = 'iL-qs (6) 

dvsv 

^qs-QR (7) 
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= qR-qp (8) 

dvpv 

^ = gp-g. (9) 

?L, Qr, Qs and gp, represent flow through the left heart, the right heart, 
the systemic capillaries and the pulmonary capillaries, respectively. The 
algebraic Equations (l)-(5) may be used to eliminate the flow variables q 
and the pressure variables p from the system of differential Equations (6)- 
(9). The result is a system of four differential equations in the four volume 
variables. 

From Assumption 1, the following is obtained: 

VSA + Vsv + VpA + Vpv = Vo (10) 

where Vq is the total blood volume, a constant. Equation (10) indicates that 
the four volume variables are not independent. As fpA has the smallest value, 
it is chosen for elimination and a system of three equations in vsa, vsv and 
vpy, is obtained (mathematically any one of the four volume variables could 
be eliminated). Hence the mathematical model of the cardiovascular system 
consists of the following system of three ordinary differential equations: 



dvsA _ Vsa ^ vsv _^ F *Cl* vpv _ vp ^^^^ 
dt Rs * CsA Rs * Csv Cpy Rs * Csv 

dvsv ^ Vsa ^ 1 ^ F*Cr ^^^^^^ ^^^^ 
dt Rs * Csa Rs * Csv Csv 



dvpv Vo - Vsa -vsv-vpv , 1 F * 

( D . ^ + -7^ PPV (13) 



dt Rp * CpA Rp * Cpv C 



PV 



Note that this model is linear in the three volume variables. The system 
becomes nonlinear, on the inclusion of the baroreflex control loop. 
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2.2.2 Model with Baroreceptor Control 
The Hill function is defined by: 

y = Ux) = (14) 

/4a;):[0,oo)^[0,l) 

where a is a constant which corresponds to the point of half the maximum 
output, /n(a) = I • As n increases the curve approaches a perfect "switch" 
or a step function at x = a. The baroreceptor response curve described in 
the literature strongly resembles a Hill function and therefore is modelled in 
this paper as: 

where B is the total baroreceptor afferent activity, n is a measure of the 
barorefiex gain, and pc, is the critical arterial pressure. The term "gain" 
normally is used to represent a ratio of the change in output to a change in 
input, for very small changes. This however is essentially the mathematical 
definition of a derivative. Thus, for our model, using the Hill function for 
the barorefiex response, we define gain by the derivative 

that is, gain yu is equal to the slope of the response function (for fixed n), 
at a particular point. For simplicity we take the value of this when psA = 
Pc, which is precisely 

J . (17) 

This is a good measure of the gain as it is very close to the maximum value 
of the slope (from calculus). 

Also, note that psA is proportional to vsa (since Csa is a constant), so B 
can be expressed in terms of vsa rather than psA- This yields: 

B„fe.)^S„(gi)^^^yi^^ (18) 
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where Vq is the volume at the critical pressure. It should be noted that 
other functions with a similar graph shape could be used for the baroreceptor 
response function, for example the hyperbolic tangent function. 

Changes in heart rate, F and systemic capillary resistance, Rs must be 
in an opposite direction to a change in arterial blood pressure, in order to 
restore normal pressure. Thus a simple model of the baroreflex action on F 
is: 

F = Foil - B)+ = (19) 

where Fq is a constant. Equation (19) implies that if B approaches 1 (i.e. 
very large pressure Psa)-, then F will be zero. However, one would expect that 
in reality F will have a non zero minimum value even when B approaches 1. 
The following is a more realistic representation: 

where Fi and F2 are constants, and F2 is the value of F when B = 1. 
Similarly, the baroreflex action on Rs is modelled as: 

R, = n.il-B)+E.= ^^r^}f)^^+Ik (21) 

On the other hand, changes in systemic venous compliance, Csv, and 
systemic venous unstressed volume, Vb, are in the same direction as a change 
in arterial blood pressure. Thus the action of the baroreflex on each of Vd 
and Csv is modelled as: 

V, = iD„B) + D. = j^^^^S^^D. (22) 

C.. = (C,.S) + C. = ^^Ji|£^ + C. (23) 

Now, from the basic linear cardiovascular model (11)(12)(13), four differ- 
ent nonlinear cardiovascular models are obtained, corresponding to insertion 
of the baroreceptor feedback function into each of F, i?s, Csv and Vd, as 
above. This allows independent investigations of each of the four feedback 
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Tabic 1: Typical parameter values for an adult human being (Hoppensteadt 
and Peskin, 1992). 



PARAMETER 


NORMAL VALUE 


Coiiipliaiice in Systeiiiic Arteries, Csa 


0.01 litres/nini Hg 


Compliance in Pulmonary Arteries, CpA 


0.00667 litres/mm Hg 


Compliance in Pulmonary Veins, Cpy 


0.08 litres/mm Hg 


Systemic resistance, Rs 


17.5 mm Hg/(litre/min.) 


Pulmonary resistance, Rp 


1.79 mm Hg/(litre/min.) 


Compliance in Right Heart, Cr 


035 litres/mm Hg 


Compliance in Left Heart, Cl 


0.014 litres/mm Hg 


Heart rate, F 


80 beats/min. 



loops, which would be very difficult to carry out in experiments on live sub- 
jects. 

2.3 PARAMETER DETERMINATION 

Many of the parameters in this model are available in the literature, as 
displayed in Table 1. The remaining parameters, Fi, F2, Ri, R2, Ci, C2, 
Di, D2, Csv, Vc and Vb are not found in the literature and need to be 
determined. 

2.3.1 Critical Volume, Vc 

No value of Vc is in the literature. However the normal resting value of vsa 
is known to be 1.0 litre. It is assumed that the resting and critical states are 
the same and hence Vc is taken as 1.0 litre. 

2.3.2 Normal Systemic Venous Unstressed Volume, Vp, 

An exact normal value of Vb is not found in the literature. However, Coleman 
(1985), gives the value of Vd as "over half" the systemic venous volume. As 
the normal value of vsy is 3.5 litres, the normal value of Vb is taken as 2.0 
litres in this study. 
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2.3.3 Normal Systemic Venous Compliance, Csv 

The normal systemic venous compliance is given as 1.75 litres/mm Hg (Hop- 
pensteadt and Peskin, 1992). This value of Csv does not account for systemic 
venous unstressed volume. However, this study considers the systemic ve- 
nous unstressed volume. Using Equation 2, with vsv — 3.5 litres, Vd = 2.0 
litres and psv = 2 mm Hg, the value of Csv is computed as 0.75 litres/mm 
Hg. 

2.3.4 Normalized Hill Function Constants, Fi, F2, Ri, R2, Ci, C2, Di, 
and D2 

These constants arc required for the use of the Hill function to model the 
baroreceptor afferent activity, B, for the systemic venous compliance, sys- 
temic venous unstretched volume, systemic resistance and heart rate. This 
current study appears to be the first time such an approach has been taken, 
which is why these constants are not available. Therefore, sets of different val- 
ues of each of these constants are investigated in this study, over ranges which 
yield baroreceptor responses consistent with experimental observations. 

3 ANALYSIS AND RESULTS 

The control of the baroreflex on heart rate, F, systemic capillary resistance, 
Rs, systemic venous unstressed volume, Vd, and systemic venous compliance, 
Csv, are investigated individually in the mathematical model. The non- 
linear baroreflex response function is substituted in each parameter in turn, 
to obtain the corresponding model for investigating the baroreflex effect on 
the parameter under consideration. See the Appendixes and [1] for details. 
We analyzed each of the four models to find out if a bifurcation occurs as 
the baroreceptor gain fi varies, using Hopf's Bifurcation Theorem. First, the 
steady state solution of each model was found. Then the system of equations 
was hnearised at this steady state. The value of the Jacobian matrix obtained 
at the steady state was found. Since it is a real 3x3 matrix, with constant 
real entries, the eigenvalues are either all three real or else one real and 
two complex conjugate. The eigenvalues of the resulting matrix were found. 
In each case there existed, for some values of /x, one real and two complex 
eigenvalues. The real part of the pair of complex conjugate eigenvalues was 
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plotted as a function of fi, to find out if a crossing point existed. The value 
of /i at which the real part crosses the /x-axis is what is known as the crossing 
point or Hopf bifurcation point. At this point there exists a pair of purely 
imaginary eigenvalues ±iuj and the steady state is said to be nonhyperbolic. 
When a crossing point was found, the imaginary part uj of the complex 
eigenvalues was plotted to obtain its value at the crossing point. The third 
(real) eigenvalue always remained negative. According to the general theory 
of Liapunov stability, when the real part of the complex eigenvalues crosses 
from negative to positive, the equilibrium state changes from asymptotically 
stable to unstable. The computation of the eigenvalues and the plotting of 
the curves was done using Maple. From the Hopf Bifurcation Theorem [8], 
generically at such a crossing point, a periodic solution is either created or 
destroyed. Further numerical computations verify the existence of a stable 
limit cycle near the crossing point. The imaginary part lv at the crossing point 
gives a good approximation to the frequency of the resulting oscillations. 

3.1 BAROREFLEX CONTROL OF HEART RATE 

Models with i?s, Csv and Vd taken as constants and F given by Equation 
20 were considered. Assuming a normal heart rate of 80 beats/min., values 
of Fi and F2 considered were: Fi = 160 beats/min. and F2 = beats/min., 
Fi = 80 beats/min. and F2 = 40 beats/min., and Fi = 40 beats/min. and 
F2 — 60 beats/min. All of these models exhibited a stable steady-state, for 
all values of gain ^ tested. No evidence of waves was found. Details of the 
models are given in [1]. 

3.2 BAROREFLEX CONTROL OF SYSTEMIC RE- 
SISTANCE 

Models with the baroreflex affecting only systemic resistance Rs, while Csv, 
Vd and F are taken as constants, are considered next. Rs is given by Equa- 
tion 21 and a typical value of systemic resistance is 17.5 mmHg/(litre/min.). 
Values of Ri and R2 used are: Ri — 35 mm Hg/(htre/min.) and R2 
— mm Hg/(litre/min.), Ri — 20 mm Hg/(litre/min.) and R2 — 7.5 
mm Hg/(litre/min.), and Ri = 15 mm Hg/(litrc/min.) and R2 = 10 mm 
Hg/(litre/min.). These models had a stable steady-state for all values of 
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gain, II tested, and showed no indications of waves. Details of the models are 
given in [1]. 

3.3 BAROREFLEX CONTROL OF VENOUS VOL- 
UME 

The baroreflex influences the systemic venous volume through the unstressed 
volume, Vd and the compliance, Csv- Models with the baroreflex control- 
ling Csv and Vb individually were considered (see Appendix for equations). 
For both cases the models exhibited unstable steady-states for gains past 
a crossing point with pure imaginary eigenvalues. Figure 1 displays four 
graphs obtained for models with the baroreflex controlling unstressed venous 
volume. It shows graphs of the real parts of the complex eigenvalues, for 
models with D2 equal to 0, 0.5, 1.0, and 1.5 litres, respectively. Note that 
Re{\) crosses through zero in all cases. This implies a Hopf bifurcation, 
giving birth to an oscillation or wave. 

Similarly, Figure 2 is obtained from models with the baroreflex control- 
ling venous compliance only. It shows graphs of the real parts of the com- 
plex eigenvalues for three cases of models with C2 equal to 0, 0.25, and 0.5 
litres/mm Hg, respectively. All three cases give a Hopf bifurcation. 

The values of the imaginary parts of the complex eigenvalues at the cross- 
ing points give the angular frequency of the oscillations produced. From these 
frequencies, the periods of oscillation of all these models were found to be 
between 7 and 12 sees. Note that this is in perfect agreement with reported 
values of Mayer waves in human subjects. 

Figure 3 shows phase portraits obtained for the system with D2 = 0, for 
gains /J equal to 10 and 20. At /i = 10, there exists no limit cycle. However, 
for // = 20, a limit cycle exists in confirmation of the above analysis. In all 
cases considered, when a crossing point existed, a limit cycle was found to 
be stable and supercritical (that is, for /i above the crossing point). 

4 DISCUSSION 

Models with only heart rate F, or systemic capillary resistance Rs, controlled 
by the baroreflex did not exhibit a Hopf Bifurcation, while the models with 
systemic venous compliance, Csv or systemic venous unstressed volume, Vo 
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controlled by the baroreflex were capable of Hopf bifurcation. Hence the 
effect of the baroreflex on F and Rs, individually, is not the cause of oscil- 
lations. However, if the effect of the baroreflex on F and Rs were combined 
in models with the effects on Vd or Csv, they may play a part in causing 
instability. 

It is observed that for all models with Vd and Csv individually controlled 
by the baroreflex, the real part of the complex eigenvalues increases as the 
gain ji increases and the graph crosses the /x-axis at a positive value of 
This implies a Hopf bifurcation and the presence of a limit cycle oscillation. 
The stability of this limit cycle oscillation has been verified numerically. The 
similarity of the results obtained for models with Csv and Vb individually 
controlled by the baroreflex is to be expected, as the two have similar effects 
on blood flow. The model with the baroreflex controlling only Csv is more 
stable than that with the baroreflex controlling only Vu. Thus it appears 
that the baroreflex control of Vd is more important than the control of Csv 
where gain induced instability is concerned, but control of Csv also plays an 
important part. Remarkably, the periods of the oscillations fall within the 
range of 7 to 12 sees, given by Penaz [18]. 

Measurements of Wesseling et. al. [21] on human subjects revealed a pe- 
riodic fluctuation in heart rate which was synchronized with the Mayer waves 
in blood pressure. In these experiments, the same period was observed for 
the modulation of the heart rate as for the blood pressure waves. The present 
model provides a possible explanation for this heart rate modulation. Assume 
that Mayer waves have arisen through the mechanism of this paper. Then 
the baroreceptors would sense a fluctuating blood pressure in the systemic 
arteries, and would therefore exert a fluctuating feedback control, affecting 
the heart rate in synchronization with the Mayer waves, as observed in [21]. 

Further insights were obtained on varying two parameters in the model 
simultaneously; namely, the gain parameter /x together with either one of 
D2 or C2 (the minimum unstressed systemic venous volume or compliance 
respectively). As either of D2 or C2, increases, the value of gain at which 
the graph crosses the /x-axis increases. This suggests that an increase in D2 
or C2 increases the stability of blood pressure. Changing the value of D2 or 
C2 causes resetting of the baroreceptor curve. This happens in an individual 
with a time constant of 10 hours and so is not of importance to this inves- 
tigation. However as different people may have different D2 and C2 values, 
different people can be expected to take different times before Mayer waves 
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are observed, when subjected to identical Mayer wave inducing stresses. In 
particular it could be expected that young and old people will have different 
D2 and C2 values and this may explain the difference in the incidence in 
Mayer waves observed in young and old people. As we would expect larger 
D2 and C2 values in older people, corresponding to veins which have become 
stretched and less fit, our observation that stability increases with increasing 
values of D2 and C2 is in agreement with the reactions observed experimen- 
tally [6, 9, 11, 14, 16]. In Figure 4, D2 is plotted against the crossing point 
value of gain /i. The top left region represents the parameter values for which 
the equilibrium state is stable, and corresponds to older subjects, who would 
tend to have smaller gains /i, and larger D2 values. The lower right region 
represents unstable equilibria, susceptible to oscillations, and corresponds to 
youth. Thus, stability depends on both the baroreflex gain fj, and D2 . Thus, 
young adults, with high gain and small D2 , arc in the unstable region, while 
older adults with the opposite characteristics are in the stable region. A sim- 
ilar situation holds for /i and C2. We conclude that the existence of Mayer 
waves and their disappearance with age can be explained, at least in part, 
as a case of gain induced instability. 
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5 APPENDIX A: CONTROL OF VENOUS 
VOLUME 

A mathematical model is constructed in which the systemic venous un- 
stressed volume Vd is controlled by the baroreflex, while Rs, Csv and F 
are assumed to remain constant. Vd is given by Equation 22. Different 
choices of the constants Di and D2 are considered. The stability of the 
equilibrium state is investigated, by computation of the eigenvalues of the 
Jacobian matrix, using Maple. 

5.1 Di = 4.0 litres and D2 = litres 

The circulation is then described by the following system of equations: 
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dvsA 4:0vsA , 8vsv , , ^ 32{vsa)'^^ 



dt 7 105 ^ "^^^ 105(1 + (^;5a)^'') 



dvsv _ ^OvsA SOvsv _^ 320{vsa) 



4m 



dt 7 21 21(1 + (^;5a)"') 



dt 



= 420 - (84^5^) - {84vsv) - {105vpv) 



The steady state values of the system are found to be: v^a = 1-0 litres, v^y 
— 3.5 litres, and vpy = 0.4 litres. Linearization of the model at the steady 
state gives the matrix A. 

/ ^40 _ 32^ _8_ 
^ - y + ^ 21 

V -84 -84 -105 

The eigenvalues of A are: 

2405 32/x 



63 315 



[/1/3 y 2405 32u I ^,,,1/, 

A2 = ^ - + -VS(U^'^ + V) 

^ 2 2 63 315 2 ^ ^ 



^ 2 2 63 315 2 ^ ^ 



where 



-5103633725 _ 3618065648/x 184832//^ ^ 32768//^ 
250047 416745 ^ 416745 ~ 31255875 ^ 

- (4488999170550 + 45082290017400//+ 



2835 

9458754238932//2 - 942218496//^ + 2248704/)^/^ 
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and 



-2876953 , 11552^ _ 1024/^^ 
Y — 3969 3969 99225 

For values of jj, of interest, eigenvalue Ai is real and negative. Eigenvalues 
A2 and A3 are complex conjugates, with real part which crosses through zero 
from negative to positive as increases, near = 18, as shown in Figure la. 
Phase portraits, on each side of the crossing point, are shown in Figure 3. 

5.2 Di = 3.0 litres and L>2 = 0.5 litres 

The circulation is then described by the following system of equations: 



dvsA 4:0vsA , Svsv , V 28{vsa)'^^ 

- + -777^ + {14:Vpv) 



dt 7 105 ' "^"^^ 105(1 + 

dvsv ^ ^OvsA SOvsv 40(1 + 7{vsa)^'') 

dt 7 21 21(1 + {vsa)'") 

dvpv 



dt 



= 420 - [Uvsa) - mvsv) - m^vpv) 



The steady state values of the system are found to be: vq,a = 1-0 litres, v^v 
= 3.5 litres, and vpy = 0.4 litres. Linearization of the model at the steady 
state gives the matrix A. 



-40 _ S^t 



14 



A= \ 4 + ^-$ 



7 ' 7 21 
-84 -84 -105 



The eigenvalues of A are: 



Ai = U^/^ -V - 



2405 
~63~ 



8yU 

105 
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U^/^ V 2405 8/x i r^.rM-^ ..x 



where 



^ _ -5103633725 904516412^ 11552^^ ^ 5^2^^ 

250047 138915 ^ 46305 ~ 1157625 ^ 

4 

-—(17955996682200 + 135246870052200;U+ 
21282197037597//2 - 1589993712/^^ + 2846016/)^/^ 



and 



-2876953 , 2888At _ 64At^ 
y _ 3969 1323 11025 



C/V3 

For values of // of interest, eigenvalue Ai is real and negative. Eigenvalues 
A2 and A3 are complex conjugates, with real part which crosses through zero 
from negative to positive as [i increases, near [i — 24, as shown in Figure lb. 

5.3 Dx = 2.0 litres and D2 = 1.0 litres 

The circulation is then described by the following system of equations: 

dvsA ^OvsA ^ 8vsv ^ , . 2A(vsa)^'' + 8 



dt 7 105 ' "^""^ 105(1 + {vsa)^'') 

dvsv ^ ^OvsA 8^vsv 80(1 + 3(^5^)'^) 
dt 7 21 21(1 + (^;5a)''') 

^ = 420 - {%Avsa) - {Uvsv) - [lOhvpv) 
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The steady state values of the system are found to be: v^a = 1-0 litres, ?;sv 
= 3.5 litres, and vpy = 0.4 litres. Linearization of the model at the steady 
state gives the matrix A. 

-40 _ 16/£ _8_ 14 \ 



A — \ 4^ 1 l^A _''8"0 n 
^ — I 7 21 21 ^ 



5 105 
Lt _8C 

84 -84 -105 ; 



The eigenvalues of A are: 

Ai = U^/^ -V- 



2405 16fi 
~63~ ~ 315 



[/1/3 Y 2405 16/i i r,rrM-^ ..x 



£V3 y_ 2405 _16^_i^ , 

^ 2 2 63 315 2^ ^ 



where 



^ _ -5103633725 _ 1809032824// 46208//^ _ 4096//^ 

250047 416745 ^ 416745 ~ 31255875 ^ 

^ -(4488999170550 + 22541145008700//+ 



and 



2835 

2364688559733//^ - 117777312/t^ + 140544//^)^/^ 



-2876953 , 5776^ _ 256^^ 
y _ 3969 3969 99225 



C/V3 

For values of /t of interest, eigenvalue Ai is real and negative. Eigenvalues 
A2 and A3 are complex conjugates, with real part which crosses through zero 
from negative to positive as /t increases, near /t = 36, as shown in Figure Ic. 
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5.4 Di = 1.0 litres and D2 = 1.5 litres 

The circulation is then described by the following system of equations: 



dvsA _ 40vsA^Svsv^.,, , 20ivsAf'^ + 12 
^ - + W + ^''"^^^ - 105(1 + (vsAr^) 



dvsv ^ A^vsA 80^gy 40(3 + 5(^5^)'^) 

dt 7 21 21(1 + (^;5a)''') 

dvpv 
dt 



= 420 - {Uvsa) - {^^vsv) - (105^;py) 



The steady state values of the system are found to be: wsa = 1-0 litres, t'sv 
= 3.5 litres, and v-py — 0.4 litres. Linearization of the model at the steady 
state gives the matrix A. 



I ' -84 -84 



7 105 105 
40 I _g0 

7 21 21 




The eigenvalues of A are: 



2405 _ 8// 
~63~ ~ 315 



\/ 2405 8^ I r^rrM-^ ..x 



where 
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^ _ -5103633725 904516412/i 11552/i2 512/^^ 

250047 416745 ^ 416745 ~ 31255875 ^ 

-47(17955996682200 + 45082290017400//+ 
2835 

2364688559733//^ - 58888656//^ + 35136//^) 

and 

-2876953 , 2888;u _ 64^^ 
Y — 3969 3969 99225 

C/V3 

For values of fj, of interest, eigenvalue Ai is real and negative. Eigenvalues 
A2 and A3 are complex conjugates, with real part which crosses through zero 
from negative to positive as fj, increases, near fj, — 71, as shown in Figure Id. 

6 APPENDIX B: CONTROL OF VENOUS 
COMPLIANCE 

A mathematical model is constructed in which the systemic venous com- 
pliance Csv is controlled by the baroreflex, while Rg, Vb and F, are held 
constant. Cgy is given by Equation 23. The stability of the model, for dif- 
ferent values of Ci and C2 is explored. The results of these calculations are 
presented in Figure 2. 

6.1 Ci = 1.5 litres/mm Hg and C2 = litres/mm Hg 

The circulation is then described by the following system of equations: 

dvsA _ ^OvsA , 4r,sr(l + (r5.0'^') , , ^ 8(1 + (r,s,0'^') 
dt - 7 + W5{vsAr' mvsA^ 

dvsv ^ ^OvsA ^Ovsyjl + (vsa)^'') 80(1 + jvsAf'') 
dt 7 21{vsa)'' 21{vsa)'' 
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^ = 420 - (8Avsa) - (SAvsv) - (105vpv) 

The steady state values of the system are found to be: vsa = 1.0 htres, 
^^sv = 3.5 htres, and vp^ — 0.4 htres. Linearization of the model at the 
steady state gives the matrix A. 

-40 _ 8fi _8_ 

7 „35 105 
40 I 80^ _80 

7 21 21 



-84 -84 




The eigenvalues of A are: 



A3 = 



63 105 



U^/'^ V 2405 8/i i 
2 2 63 105 2 ^ ^ 



where 



-5103633725 _ 904516412// 11552//^ _ 512//^ 
250047 138915 ^ 46305 ~ 1157625 ^ 

4 

;i7955996682200 + 135246870052200//+ 



and 



2835 

21282197037597/i2 - 1589993712//3 + 2846016//^)^/^ 



-2876953 , 2888jU _ 64;u^ 
Y _ 3969 1323 11025 



C/V3 

For values of /t of interest, eigenvalue Ai is real and negative. Eigenvalues 
A2 and A3 are complex conjugates, with real part which crosses through zero 
from negative to positive as /i increases, near n = 24, as shown in Figure 2a. 
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6.2 Ci = 1.0 litres/mm Hg and C2 = 0.25 litres/mm 
Hg 

The circulation is then described by the following system of equations: 



dvsA ^OvsA , Svsv{1 + {vsa)''') , , 16(1 + 



dt 7 35(1 + 5{vsA)'n ' ' 35(1 + 5(^5a)'^) 

dvsv ^ /^QVSA ^Qvsv{1 + {vsaY^) + {vsAf^) 

dt 7 7(1 + 5(^;5a)'^) 7(1 + 5(^5a)"^) 

dvpv 



dt 



420 - {Uvsa) - {Uvsv) - (105^;pv) 



The steady state values of the system are found to be: vsa = 1-0 litres, 
vsv — 3.5 litres, and t'pv = 0.4 litres. Linearization of the model at the 
steady state gives the matrix A. 

_ 16i£ _8_ -14 \ 

7 105 105 
40 I 16D/i 80 Q 



7 ' 21 21 

V -84 -84 -105/ 



The eigenvalues of A are: 



63 315 



C/V3 V 2405 16/i I r,rrAn t.x 

\ ^ + - V3([/^/^ + \/) 

2 2 63 315 2 ^ ' 



K_ 2405 _16,_. 
^ 2 2 63 315 2^ ' 



where 
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^ _ -5103633725 _ 1809032824^ 46208^^ ^ 4096/i3 
~ 250047 416745 ^ 416745 ~ 31255875 ^ 

-(4488999170550 + 22541145008700^^+ 



2835 

2364688559733//^ - 117777312//=^ + 140544//^)^/^ 

and 

-2876953 , 5776/f _ 

3969 3969 99225 



V 



For values of /x of interest, eigenvalue Ai is real and negative. Eigenvalues 
A2 and A3 are complex conjugates, with real part which crosses through zero 
from negative to positive as ji increases, near = 36, as shown in Figure 2b. 

6.3 Ci = 0.5 litres/mm Hg and C2 = 0.5 litres/mm Hg 
and 

The circulation is then described by the following system of equations: 

dvsA^ ^QVSA Avsv{l + {vsAf^) , 8(1 + [vsAf^) 

dt 7 ^ 35(l + 2(^^^)^'^) '^''^^ 35(1 + 2(^^^)^0 

dvsv ^OvsA ^Ovsvil + (vsa)^^) ^ 80(1 + Vsa^^ 



dt 7 7{1 + 2{vsa)'^) 7{l + 2{vsA)'n 

dvpv 



dt 



420 - {Uvsa) - {^^vsv) - {l^hvpv) 



The steady state values of the system are found to be: v^a = 1-0 litres, 
t'sv = 3.5 litres, and v-p\ = 0.4 litres. Linearization of the model at the 
steady state gives the matrix A. 

/ ^40 _ 8iJ_ 8 



7 105 105 
40 I 80^ _80 

7 21 21 



v 
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The eigenvalues of A are: 



A2 = 



63 315 



y 2405 8u i r,rrAn ..x 



^ 2 2 63 315 2 ^ ^ 



where 



-5103633725 904516112// 11552//2 512/r 



250047 416745 416745 31255875 

4 , 

17955996682200 + 45082290017400//+ 



and 



2835 

2364688559733//^ - 58888656//^ + 35136//^) 



-2876953 , 2888^* _ Mj? 
Y _ 3969 3969 99225 



C/V3 

For values of /t of interest, eigenvalue Ai is real and negative. Eigenvalues 
A2 and A3 are complex conjugates, with real part which crosses through zero 
from negative to positive as /t increases, near fj, — 71, as shown in Figure 2c. 
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